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ABSTRACT 

We obtained Spitzer/TRAC 3.6 - 8 /jm images of the nearby spiral galaxy NGC 4258 to study 
possible interactions between dust and the radio jet. In our analysis we also included high-resolution 
radio continuum, Ha, CO, and X-ray data. Our data reveal that the 8 /Ltm emission, believed to 
originate largely from PAH molecules and hot dust, is an excellent tracer of the normal spiral structure 
in NGC 4258, and hence it originates from the galactic plane. We investigated the possibility of 
dust destruction by the radio jet by calculating correlation coefficients between the 8 ^m and radio 
continuum emissions along the jet in two independent ways, namely (i) from wavelet-transformed 
maps of the original images at different spatial scales, and (ii) from one-dimensional intensity cuts 
perpendicular to the projected path of the radio jet on the sky. No definitive sign of a correlation (or 
anticorrelation) was detected on relevant spatial scales with either approach, implying that any dust 
destruction must take place at spatial scales that are not resolved by our observations. 
Subject headings: galaxies: jets — galaxies: ISM — galaxies: individual (NGC 4258) — galaxies: 
spiral — infrared: galaxies — radio continuum: galaxies 



1. INTRODUCTION 

The nearby (D=7.2 Mpc; iHerrnstein et all 11999ft 
spiral galaxy NGC 4258 is known for its cen- 
tral black hole ([Mivoshi et al.l 119951 : IHerrnstein et al.l 
11998ft and jets (anomalous arms), and is unique 
in its geometry among nearby spiral galaxies. It 
possesses kpc-scale anomalous radio continuum and 
Ha and X-ray arms, in projection perpendicu l ar to 
its normal spiral arms (ICourtes fc Cruvellierl HMD: 



van der Kruit. Port, fe Mathewsonlll972ft~ that have been 
shown to connect to a nuclear jet with sub-arcsecond res- 
olutio n observations (jHerrnstein et al.l 119971 : iCecil et all 
2000). The Ha em ission along the jet is no t ex- 
cited by OB stars (jCourtes fe Cruvellierl 11961ft but 
instead is consistent with shock excitation (e.g., 
ICecil. Morse, fe Veilleuxl 11995ft . The anomalous radio 
arms emerge as linear features from the nucleus up to a 
projected distance of 840 pc, and then have a first bend 
of about 40°, followed by further bends , in the trailing 
sense relative to the galactic rotation ()Krause fe Lohrl 
I20i51 . 

The accretion disk of the central engine has an incli- 
nation angle of 83° and a position angle (p. a.) of 86° 
(jMivoshi et al.lll995l ). i.e., it is oriented nearly east-west. 
The p. a. of the galactic disk, however, is 15 0°, which is 
nearly north-south, and its inclination is 72° (|van Albadal 



1980). The plane of the galactic disk and the plane 
of the accretion disk have significantly differing orienta- 
tions, being tilted by 60° - 83° with respect to each other 
(jKrause. Fendt. fe Neini nger 200.3). As the jets emerge 
perpendicular to the accretion disk, they have to travel 
through the galactic disk at a rather small angle. The 
distribution of the molecular clouds (see below) and the 
Ha emission of the jets indicate a direct interaction of 
the jets with disk material within 2 kpc on either side of 
the nucleus. 

Several studie s (e.g., ICecil. Wilson, fe Tullvl U992I: 
ICecil et aT1l2000l ) have elaborated on the "jet" structure 
of the anomalous radio continuum arms. One of the most 
comprehensive studies in the radio contin u um a nd opti- 
cal wavelengths was made bv lCecil et al.1 (pOOOft . They 
find the jet to be in projection perpendicular to a nu- 
clear accretion disk (but in space misaligned by 65°), 
and associated with two radio hot spots at around 1-2 
kpc from the nucleus. They find optical emission-line 
arcs on the leading edges of the radio hot spots. They 
explain the observed structure with a shock that prop- 
agates with a velocity of 350±100 km s _1 but is now 
stalled. The shock connects to the anomalous radio con- 
tinuum arms. They interpret the overall structure with a 
precessing jet that has recently moved out of the plane of 
the galaxy. The outer anomalous arms may represent gas 



2 



Laine et al. 



that has been driven out of the plane by the jet-affected 
halo gas that is now impacting the disk, causing an ap- 
parent curvature when projec ted to the plane of the sky 
(IWilson. Yang, fc Ceci]ll200l . 

Single-dish 12 CO 1-0 line observations carried out 
with the 30-m IRAM telescope revealed molecular gas 
in NGC 4258 only in t he center and along the Ha jets, at 
dista nces up to 2 kpc (jKrause et al.lll990UCox &; Downesl 
1990). It was not detected in other parts of the (cen- 
tral) galactic disk, down to the detection limit of about 
30 mK brightness temperature. The velocities of the 
molecular gas correspond primarily to the rotation of 
the galactic disk. Hence the molecular clouds are in 
the plane of the disk. Hi g h-reso lution 12 CO 1-0 ob- 
servations of iKrause et al.l (|2007j ) with the IRAM in- 
terferometer at Plateau de Bure revealed that the bulk 
of the CO emission comes from two parallel ridges on 
both sides of the linear radio jets that cross the nu- 
cleus. This morphology was confirmed by BIMA 12 CO 1- 
and Sub-Millimeter Arra y (SM A ) 12 CO 2-1 q bserva- 
tions (|Sawada-Satoh et al.ll2007t ). Krause et all (|2007l ) 
proposed a model in which the molecular gas temporar- 
ily accumulates along ridges surrounding the jet due to 
an interaction of the jet magnetic field with the molecular 
gas clouds. 

While star formation has presumably not (yet) been 
triggered by the radio jet, it is of interest to study how 
the dust and gas have been affected by the burrowing 
radio jet. Previou s studies of galactic jet-g as or jet-dust 
interactions (e.g.. iWhittle fc Wilson! [20041) suggest that 
molecular gas clouds that cross the path of the jet will 
cause bends and knots in the jet, but eventually the jet 
may penetrate and ablate the cloud, in the process ioniz- 
ing and accelerating the gas. In radio galaxies the high- 
velocity jets (fjet = 10 4 - 10 5 km s _1 ) coming from the 
nuclei are expec ted to destroy dust grains very quickly 
(IDe Youndll998h . In NGC 4258 the jet speed is an order 
of ma gnitude smaller, around 1000 km s" 1 (Cec il et al.l 
1992), and therefore it may be insufficient to cause dust 
destruction. 

We obtained Spitzer /Inir&red Array Camera (IRAC) 
images of NGC 4258 from 3.6 to 8 /jm. The 8 fim image 
is dominated by 7.7 fim and 8.6 fim polycyclic aromatic 
hydrocarbon ( PAH) emission feat ures and radiation from 
hot dust (e.g. JRegan et af1l2004[ ). Similarly, the 5.8 fim 
image is dominated by the 6.2 /im PAH emission feature. 
We compare the IRAC images to images of the radio 
continuum, CO, Ha, and X-ray emission at a compara- 
ble resolution to investigate how the radio jet may have 
affected the dust properties. We use the 8 fim emission 
as a tracer of dust emission, the 3.6 cm radio continuum 
and X-ray emission as a tracer of the jet, and we look for 
signs of possible dust destruction by the jet by search- 
ing for an anticorrelation between these emissions. We 
do this with the help of a wavelet analysis and by look- 
ing at one-dimensional intensity slices taken across the 
projected position of the jet in the images. 

2. OBSERVATIONS 

NGC 4258 was observed with IRAC dFazio et al.H200l 
aboard the Spitzer Space Telescope on 2005 December 24. 
Observations with the channel 2 (4.5 /im) and 4 (8 fim) 
cameras were centered on the galaxy. Observations in 
channels 1 (3.6 fim) and 3 (5.8 /ion) also covered the 



central part of the galaxy. The total mapped area was 
about 28' x 26' in each channel. We used the small- 
scale cycling dithering option together with half array 
map steps in the array column and row directions to 
cover the required area on the sky and to produce suf- 
ficient redundancy against data artifacts. We also used 
the high-dynamic-range observing mode which takes a 
short-exposure image before a long-exposure image. The 
short-exposure frame times were 1.2 seconds and the 
long-exposure frame times were 30 seconds. In this pa- 
per we use only the 30-second frames, as we did not find 
any saturated areas in the galaxy. The pixel size in each 
channel is roughly l'/2 x l'/2. The observations reported 
here use the S14 pipeline processing version. We ran 
the individual IRAC basic calibrated data (BCD) frames 
through the IRAC artifact mitigation software (available 
as contributed software on the Spitzer website). Subse- 
quently we improved the muxbleed artifact mitigation in 
the long-exposure images from channels 1 and 2 by cus- 
tom software, kindly provided by J. Hora. We then made 
mosaics of the long-exposure frames with the Mosaicking 
and P oint Source Extraction ( MOPEX) mosaicking soft- 
ware fMakovo z fc Khanl 120051 ). available at the Spitzer 
Science Center website. We corrected for the background 
variations in the individual frames by running the BCDs 
through the overlap correction algorithm first. This al- 
gorithm normalizes the spatially overlapping parts of in- 
dividual BCD frames to the same mean value. The final 
pixel size in the mosaics is 0'.'6 x 0'.'6. MOPEX rejects 
outliers, such as radiation hits, during the construction 
of the mosaic. The standard interpolation method was 
used, whereby the relative overlapping areas of the input 
pixels from the individual BCDs are used as weights for 
calculating the value of the output pixel. 

3. RESULTS 
3.1. Morphology 

The final IRAC mosaics are shown in Figure [1] The 
emission at 3.6 fim and 4.5 fim is very smooth over most 
of the galaxy, and traces the older stellar population in 
the disk. The strongest star formation regions in the spi- 
ral arms stand out from the smooth overall distribution 
of the emission at these wavelengths. There is likely to be 
a contribution from bright supergiant stars and possibly 
from hot dust in regions of strong recent or active star 
formation. The active nucleus of NGC 4258 shows up 
as a bright, point-like component, which is likely to have 
a larger contribution from hot dust as the wavelength 
increases from 3.6 fim to 8 /im. The morphology be- 
comes patchier in the 5.8 and 8 /tm images. At 5.8 /im, 
the stellar contribution decreases dramatically and the 
6.2 /an PAH and hot dust continuum components be- 
come stronger. The 8 fim IRAC band still traces hot 
dust, but the PAH emission is likely to dominate. The 
hot dust and PAH emissions appear to trace the spiral 
arms very well (as discussed in the next section) and are 
found around star formation regions, but the anomalous 
Ha and radio continuum arms have very little hot dust 
or PAH emission. A similar change in morphology as the 
wavelength increases from 3.6 /im to 8 u m has been seen 
before, e.g., in M81 (|WiUner et al.|[200l . 

We also measured the ellipticity in the 3.6 fim and 4.5 
fim images at radii within the major star forming spiral 
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Figure 1. Grayscale 3.6, 4.5, 5.8, and 8 (im IRAC images of NGC 4258. 



arms. The ellipticity rises to a value of around 0.35-0.4 
at semimajor axes from ~ 3" to ~ 30" (100 pc - 1 kpc). 
The position angle is around 146° - 148°. This is very 
similar to the 150° position angle and 72° inclination of 
the galactic dis k of NGC 4258, as determined from H I 
observations by Ivan Albadal (fl980h . We do not see any 
evidence for a lense-like structure (bar) at a position an- 
gle of ~ 10° (contrary to what was found by Ivan Albadal 
119801 in H I) in our images, but we cannot rule out the 
existence of an oval distortion either. 

We do not see the radio continuum hot spots (about 50 
arcseconds north and 25 arcseconds south of the nucleus) 
in the 8 jum emission. They are located in relative low 
intensity valleys of the 8 micron emission. The predicted 
flux densities of the hot spots, as extrapolated from the 



radio continuum and assuming a spectral index of 0.7, 
would be less than 1 /xjy, well below the extended dust 
continuum background radiation from the disk which has 
fluxes about 2.5 mJy at the locations of the radio hot 
spots (as integrated in 5" x 5" boxes around the hot 
spots). 

3.2. 8-micron Emission as a Tracer of Spiral Structure 

We compared the maps at 5.8 /im and 8 jum with vari- 
ous other high-resolution observations of NGC 4258 that 
were available to us, such as the 3.6-cm radio continuum 
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Figure 2. Grayscale 8 /mi im age of NGC 4 2 58 wi th overlaid con- 
tours of the He? emission from lKrause et al.l (f2007l '>. The contours 
are in arbitray units and logarithmically spaced at 10.0 X (2, 4, 8, 
16, 32, 64, 128). 
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Figure 3. Grayscale 8 fim image of NGC 4258 with overlaid con- 
tours of the high resolution radio continuum em ission at 8.44 GHz 
(HPBW = 2'/2 X 2'/4) from lKrause fc Lohril200l The contours are 
at 6.0 X HT 6 X (4, 6, 8, 16, 32, 64) Jy beam" 1 . 
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Figure 4. CO emission shown with contours on a grayscale image 
of the 8 fim emission. The contours are shown at 2.0 X (3, 4, 6, 8, 
10, 12, 14, 16, 18, 20, 22, 24) Jy beam -1 km s _1 . The radio jet 
position angle is —3° for r < 24", and —43° for r > 24", whereas 
the CO and Ha (and X-ray) emission are aligned at about —25° . 

radio continuum in Figures [5] and [31 respectively. As 
Figure [2] shows, there is a very good spatial correspon- 
dence between the 8 /zm and Ha emission regions along 
the normal spiral arms of NGC 4258. However, there is 
much less spatial correspondence in the central part of 
the galaxy and along the jets. The spatial correspon- 
dence of the 8 /Ltm emission with the 3.6-cm radio con- 
tinuum emission (Fig. [3j is much weaker. They only 
coincide around the nucleus and along the star-forming 
areas of the normal spiral arms. It is known that the ra- 
dio continuum emission outside the normal star forming 
spiral arms is mainly synchrotron emis sion from the jets 
(e.g.. iHummel. Krause. fc L esch 1989;), whereas the Ha 
emission t races both the no rmal spiral structure and the 
jets (e.g.. iCecil et~aT] [l992). The X -ray emission mor- 
phology is similar to the radio continuum morphology, 
although in the central two arcminutes the jet-like fea- 
ture in the X-ray image is almost pe r fectly straight, as 
can be seen in Fig. 1 of lYang et al.l (|2007t ). At larger 
distances the bulk of the x-ray emission follows the radio 
continuum very closely, although the relative intensities 
at various positions differ. The Ha emission along the jet 
follows the X-ray em ission closely. It has been suggested 
(jWilson et al.ll2001[ ) that the jet is inclined by about 30° 
to the disk of the galaxy and exits the disk at a distance 
of a few hundred pc from the nucleus. Beyond the nu- 
cleus the linear section of the X-ray jet could be due to 
shocks, triggered by the out-of-the-plane radio contin- 
uum jet, impacting the disk plane and heating the gas 
there. The curvature of the X-ray and radio continuum 
jets beyond a few kpc from the nucleus could then be 
due to the lower density shock-heated gas that has been 
pushed out of the galactic plane (to the other side of the 
galactic plane with respect to the ra dio jet), thus causin g 
the curvature, as seen in projection (| Wilson et al.ll200lh . 

Quite a different spatial distribution is visible in 
12 CO 1-0 line observations when compared to the 8 /j,m 
emission (Figure @| . The CO emission is aligned along 
the jets and concentrated within the innermost 2 kpc in 
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radius, unlike what is seen in the 8 /tni dust emission. 
The alignment of the CO emission on the sides of the jet 
was interpreted as the molecular clouds of the disk inter- 
acting with the magnetic field of the jets (jKrause et al.1 
2007). Such an interaction is, however, improbable for 
hot dust or PAHs. 

We conclude that the 8 /im emission, emitted mainly 
by PAH dust features and hot dust (which are evidently 
components of the galactic disk), is an excellent tracer 
of the normal spiral structure of NGC 4258. This is sup- 
ported by the wavelet analysis of the various emissions, 
described in the following section. 

3.3. Wavelet Analysis of the 8-micron, Radio 
Continuum, And Ha Components in NGC 4258 

To investigate the origin of the 8 /xm emission and its 
connection to the radio and Ha emissions, we compared 
the emitting structures on different spatial scales using 
a wavelet transformation. This allows us to detect sys- 
tematic correlations or anti-correlations on each spatial 
scale. 

The wavelet technique allows us to analyze the power 
emitted at different spatial scales (through the wavelet 
spectru m) , and is more r obust against noise than Fourier 
filtering (|Frick et al.ll200ll ). Cross-correlating the wavelet 
spectra from each wavelength lets us compare the emis- 
sion structures in the various images systematically as a 
function of the spatial scale. We apply a two-dimensional 
wavelet transformation to separate the diffuse emission 
components from compact structures in the 8 /im, 3.6 cm 
radio continuum, and Ha images. The wavelet coef- 
ficients of a 2D continuous wavelet transformation are 
given by 

i r + °° x' - x 

W(a,x) = — /(x')V*( )rfx', (1) 

where ip(x) is the analyzing wavelet, x = {x,y) posi- 
tional coordinates, /(x) is a two-dimensional function 
(the image being analyzed) , and a and k are the scale and 
normalization parameters, respectively, (th e *-symbol 
denot es the complex conjuga t e). F ollowing iFrick et all 
(120011) and iTabatabaei et all (|2007[ ). we use the "Pet- 
Hat" wavelet to obtain a good separation of scales and 
to find the scale of the dominant structures in the galaxy. 
The "Pet-Hat" wavelet is defined in Fourier space by the 
formula: 

m = { cos2 (f lo & (&)) * < k < ^ (2 ) 

[ k < it or k > 47r, 

where k is the wave vector and k = |k|. The above trans- 
formation decomposes an image into "maps" of different 
scale. In each map, structures with the chosen spatial 
scale are prominent as they have larger coefficients than 
those with smaller or larger spatial scales. 

After converting the 8 /im, 3.6 cm, and Ha images into 
the same geometry, size, and resolution (2'.'5), we de- 
composed them into 10 spatial scales ranging from 5" to 
165'Q using Equations [T] and [5] Figure [5] shows the orig- 

1 To have both physically meaningful results and a sufficiently 
large number of independent points, the scale a varies between a 



inal 8 /im, 3.6 cm, and Ha images together with their de- 
composed maps on three scales. On the scale of a = 7.7" 
(270 pc), the 8 /im emission reveals small clumps that 
make up segments of rings or spiral arms, even where 
the Ha emission is faint or absent. Furthermore, in con- 
trast to the Ha and radio continuum emissions, there 
are filamentary and flocculent structures at 8 /im, some 
of which appear more extended in the northern than in 
the southern half of the galaxy. The outer spiral arms 
are better visible on the a = 16.5" (580 pc) scale. At 
this scale the strong star-forming region in the northern 
arm is the only common structure with the radio emis- 
sion. On the next larger scale, a — 35.5" (1.2 kpc), the 
dominant dust emitting structures are found in the cen- 
ter and around the strong star-forming region located in 
the northern spiral arm. 

Apart from the center the radio continuum emission is 
dominated by the jet and the star-forming region in the 
northern spiral arm on all scales. A large amount of the 
Ha emission appears as dispersed clumps of star-forming 
regions. The Ha emission jet appears as a smooth, elon- 
gated structure on scales larger than 0.5 kpc. However, 
there is not any trace of a jet-like structure in the 8 /im 
wavelet-decomposed maps. The 8 /im emission traces the 
stellar spiral arms much better than the Ha or radio con- 
tinuum emissions. 

To quantify any correlation or anti-correlation of the 
dust emission with the Ha and radio continuum emis- 
sions along the jet, we focus on the region covered by 
the radio jet. We made maps of the radio continuum, 
Ha, and 8 /mi emissions by including only regions corre- 
sponding to the radio continuum intensities larger than 
a threshold (the point sources and the stellar spiral arms 
had been subtracted from the radio continuum map be- 
fore doing this). These maps were then decomposed in 
the same way for the whole galaxy. In Figure [6] we plot 
the wavelet spectrum, i.e., the distribution of emitting 
power in terms of the spatial scale for each map given by 

/ + oo r + oo 
/ \W(a,x)\ 2 dx. (3) 
-oo J — oo 

A smooth distribution of the power spectrum is ex- 
pected if the emission along the masked radio jet orig- 
inates from a real jet. The radio emitting power in- 
creases toward the 20" - 40" (0.7 - 1.4 kpc) scales and 
decreases smoothly beyond these scales. There does not 
appear to be any overall dominant power at small spatial 
scales. However, the smallest analyzed spatial scale (5" 
or 170 pc) has the most power among the first four ana- 
lyzed spatial scales at 8 /tm. This indicates patchiness of 
the dust emission along the radio jet region. These emis- 
sion patches probably emerge from gas and dust clouds 
in the disk. It is notable that all three wavelengths share 
the same broad peak in their wavelet power spectrum 
between 20" and 40". This reflects mostly the enhanced 
emission (and size of the) central region. 

The wavelet spectrum of the Ha emission is intermedi- 
ate in nature between those of the radio and dust emis- 
sions. It is not dominated by the small scale emission, 
in contrast to the 8 /tm emission, but it also does not in- 
minimum of about twice the resolution and a maximum of about 
half of the image size. 
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Figure 5. From top to bottom: IRAC 8 fim, Ha, and 3.6cm radio continuum maps at 2'.'5 resolution (first column) and their wavelet- 
decomposed maps on a scale of 7" 7 (second column), 16" (third column), and 35'.'5 (fourth column). At the distance to NGC 4258, 1" 
~ 35 pc. The maps are shown in the north up, east left orientation and are centered on the nucleus of the galaxy. 
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Figure 6. Wavelet power spectrum for the 8 ^m, Ha, and 3.6 cm 
emissions along the track covered by the radio continuum jet. The 
data points correspond to the scales of 5.2, 7.7, 11.3, 16.6, 24.3, 
35.6, 52.3, 67.7, 112.5, and 165.0 arcseconds. The power is shown 
in arbitrary units. 

crease smoothly toward the maxima at 0.7-1.4kpc scales 
(unlike the radio continuum power spectrum). Ha shares 
the smooth decline of the other two power spectra past 
the peak. This may be a consequence of the Ha emis- 
sion consisting of both emission from the Ha jet and the 
underlying disk structures, specifically H II regions. 

We calculate the correlation coefficients on each scale 
with 

_ //^(o.x) Wr( a ,x)Ac 

where the subscripts refer to two images of the same size 
and linear resolution. The value of r w varies between — 1 
(total anticorrelation) and +1 (total correlation). 
The uncertainty in r w can be approximated by 



Ar w (a) 



(5) 



Figure 7. Wavelet cross-correlation of the 8 (im emission with the 
Ha and 3.6 cm radio continuum emissions along the track covered 
by the radio continuum jet, as a function of spatial scale. 

where n = 2.13(L/a) 2 is the number of independent 
points and L is the size of the maps. Plotting r w against 
the spatial scale shows how well structures at differ- 
ent scales are correlated in intensity and location. The 
wavelet cross-correlation coefficients r w are shown in Fig- 
ure [7J The dust emission is not correlated with the radio 
emission at small scales. A significant 8 /mi - 3.6 cm cor- 
relation (r w > 0.750 although with large uncertainty), 
is found at large scales (a > 36") due to the extended 
emission from the bright central part and because of the 
artificial smearing of the two emissions by using the same 
masking (at 2a radio continuum rms level, covering the 
jet) in both maps. In general, the dust emission has a 
better correlation with the Ha emission (even in the ra- 
dio jet region) than with the radio continuum emission. 
The weakest correlation occurs around a at about 

2 A coefficient r w of > 0.75 means that more than 50% of the 
total variation in one of the variables can be explained by a linear 
relationship between the two variables. 
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the width of the Ha jet. A good small spatial scale 8 /zm 
- Ha correlation is possibly due to the star-forming re- 
gions in the disk that are projected on the radio-jet along 
the line of sight. If so, then a superposition of the emis- 
sion from the jet and the underlying galaxy disk makes 
it very difficult to detect an anticorrelation between the 
dust and the radio emissions along the jet. 

The wavelet analysis confirms that the 8 /im emission 
traces the normal stellar spiral arms much better than 
the radio continuum and Ha emissions. The long con- 
tinuous spiral arm in the northeast stands out, as does 
the split stellar spiral arm in the south and a long section 
of a spiral arm in the northwest. It is also clear that the 
emission just north (northeast) of the nucleus at declina- 
tion 47° 19' comes from a spiral arm rather than from the 
jet. In this region no CO line emission has been detected, 
which is consistent with the result t hat the CO emissio n 
is only found flanking the jet (see IKrause et al.l fi990). 
Even more notably, the wavelet analysis shows that the 
8 /jm emission does not trace the radio continuum jet 
emission at all. 

3.4. Emission Profiles Across the Jet 

In this section we further examine the dust to radio 
continuum correlation by taking one-dimensional slices 
in the 8 /im and 3.6-cm radio continuum emissions per- 
pendicular to the sky-projected orientation of the radio 
jet at several positions across the galaxy. The positions 
of the slices were chosen to sample the radio jet at various 
positions along its length. The slice positions are shown 
in Figure [51 and the emission profiles along the slices in 
Figure [5J The only slice in which visual inspection shows 
evidence of an anticorrelation between the radio contin- 
uum and 8 /iin emissions is at position 10, which lies far 
out in the outer spiral arms, where the radio continuum 
jet has split into several branches. We believe that this 
is the result of a fortuitous distribution of the two emis- 
sions at this location, rather than of any physical PAH 
grain destruction by the radio jet. In several other slices 
we see a phase offset between the radio continuum and 8 
/im emissions, but there does not appear to be any clear 
pattern in the magnitude of the offset, in the leading or 
trailing sense, as a function of the distance from the nu- 
cleus of the galaxy or the side of the galaxy (and jet). 
We varied the orientation and the length of the slices in 
a few cases and confirmed that the lack of any correlation 
or anticorrelation in the slices is not due to the orienta- 
tion or length of the slices (see also the note below about 
performing the correlation coefficient calculation by ex- 
cluding points at the ends of the slices). 

We can quantify the preceding remarks, obtained from 
visual inspection, by calculating correlation coefficients 
between the 8 /im and the radio continuum emission 
profiles for each slice presented in Fig. |H1 We calcu- 
lated the parametric Pearson (linear) correlation and 
the non-parametric Spearman (rank-order) correlation, 
along with their P- values, for all the slices. In princi- 
ple, the Spearman correlation is more appropriate since 
there is no reason to expect that the emission profiles 
are drawn from normally distributed populations. In 
practice, both methods produced similar results, which 
are shown in Table [TJ The P- value is the probability of 
obtaining the calculated correlation coefficient assuming 
the null hypothesis, i.e., that the correlation coefficient 



is in fact zero. If this probability is lower than a con- 
ventional value, often taken to be 0.05, the correlation 
coefficient is considered to be statistically significant. 

Using a threshold of ±0.7 for a correlation or anticor- 
relation, respectively, we find only 2-3 slice profiles out 
of 15 that fulfill this criterion. Changing the threshold to 
±0.6 does not increase the number of possible correlated 
or anticorrelated slice profiles by more than 1. Slice 4 
has the most correlated profiles, and they are positively 
correlated. The next most correlated slices are 1 (nega- 
tively) and 2 (positively) , followed by slices 5 (positively) 
and 10 (negatively). There is one more positively than 
negatively correlated slice, but this is not a statistically 
significant result. 

We note that it is possible to obtain higher correlation 
numbers by excluding points at the ends of the slices. 
However, when doing so, we still do not find strong evi- 
dence of an anticorrelation (or correlation) anywhere else 
except at slice position 10. 

We repeated this correlation analysis for the 8 /im and 
X-ray images and found no evidence of an anticorrela- 
tion (or correlation) between those two emissions. The 
only position where we see a significant anticorrelation 
is a little northwest of slice position 9 in Figure [5J where 
the radio continuum emission is relatively weak (but X- 
ray emission has a stronger ridge). However, this "anti- 
correlation" is due to the one-sided distribution (on the 
northeastern side only) of the 8 /im emission with respect 
to the X-ray jet, and therefore does not provide evidence 
for any destruction of PAH molecules and dust in the jet 
path. 

4. DISCUSSION AND CONCLUSIONS 

Our data reveal that the 8 /im emission (and the 5.8 
/im emission), believed to originate largely from PAH 
dust features and hot dust, is an excellent tracer of the 
normal spiral structure in NGC 4258, and hence it origi- 
nates from the galactic plane. It allows us to distinguish 
between the normal spiral arms and the jet, which is 
not trivial, especially in the northern part. The other 
well-known spiral arm tracer, the Ha emission, cannot 
be used to distinguish the spiral arm and jet emissions 
in NGC 4258 as we find Ha emission (shock ionized) also 
along the inner jets. 

We searched for a possible signature of dust destruc- 
tion by the large-scale jets by comparing the 8-micron 
emission to radio continuum emission in two different 
ways. A wavelet analysis shows that the 8 /im emission 
exhibits a patchy structure, even along the large-scale 
jet, unlike the emission seen in the radio continuum (or 
Ha). All three studied emissions (8 /im, Ha, and radio 
continuum) show largest wavelet powers at around 30" 
(~ 1 kpc) scales, mostly due to the central region. A 
correlation analysis of the wavelength-transformed emis- 
sion at various scales shows that the 8 /im emission is 
not correlated (or anticorrelated) with the radio contin- 
uum emission on any scales. Overall, the wavelet analysis 
shows that the 8 /im emission traces the normal spiral 
arms much better than the radio continuum or Ha emis- 
sions. 

We also took one-dimensional intensity profile cuts per- 
pendicular to the radio continuum jet at several arbitrary 
positions along the jet. We only detected an anticorre- 
lation between the 8 /im and 3.6-cm radio continuum 
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Figure 8. Same as Fig.[3]but with the labeled positions of the slices taken across the radio jet and shown in Figure l9l 



emissions far out in the disk where the jet has several 
branches. We consider this anticorrelation to be a for- 
tuitous superposition of the two emissions rather than 
evidence for dust destruction. The profiles at smaller 
radii along the jet sometimes show a phase shift between 
the two emissions, but no systematic offset was found as 
a function of the distance from the nucleus or the side of 
the jet (north or south) with respect to the nucleus. We 
quantified these results by a thorough correlation coeffi- 
cient analysis. Similar results (no correlation or anticor- 
relation) were obtained when comparing the 8 /im and 
X-ray emissions. 
In conclusion, we found that the large-scale jet speed in 



NGC 4258 (likely around 1000 km s-^ lCecil et"aHH99l 
is not sufficiently high to destroy dust traced by PAH 
emission on scales which are resolved by our observations. 
Future observations with higher resolution may detect 
dust destruction on smaller scales. 
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Figure 9. Slices taken across the radio jet as shown in Fig. [Sj in 8 fim (solid line) and 3.6 cm radio emission (dashed line). 
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Table 1 

Correlation coefficients between the 8 /^m and radio 
continuum emissions along the slices shown in Figure [9] 



Pearson Spearman 
Slice Correlation Correlation 





Coefficient 


P-value 


Coefficient 


P-value 


1 


-0.68 


4.3 X 10" 5 


-0.75 


2.3 x 10" 6 


2 


0.88 


1.2 X 10~ 12 


0.75 


1.1 x 10" 7 


3 


0.41 


1.5 x 10" 2 


0.49 


2.8 x 10~ 3 


1 


0.89 


1.7 x 10~ 13 


0.91 


9.9 x 10 -15 


5 


0.55 


3.8 x 10" 4 


0.61 


7.2 x 10~ 5 


6 


-0.22 


1.8 x lO" 1 


-0.27 


1.0 x 10" 1 


7 


0.22 


2.0 x lO" 1 


0.25 


1.4 x lO" 1 


8 


0.25 


6.8 x 10 -2 


0.41 


2.5 x 10 -3 


9 


0.16 


3.0 x 10- 1 


0.08 


6.1 x lO" 1 


10 


-0.58 


7.4 x 10~ 4 


-0.60 


4.2 x 10~ 4 


11 


0.17 


1.9 x 10" 1 


0.22 


9.0 x 10~ 2 


12 


0.32 


4.0 x lO- 3 


0.35 


1.8 x 10" 3 


13 


-0.27 


3.7 x 10" 2 


-0.08 


5.4 x 10" 1 


14 


-0.32 


1.8 x 10" 2 


-0.44 


6.8 x 10" 4 


15 


0.30 


1.0 x 10" 1 


0.35 


5.4 x 10" 2 
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